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^This  paper  examines  the  cyclotron  marrer  1 it 


general  magnetic  harmonic  number  transverse  electric  (TE) 
and  transverse  magnetic  (TM)  waveguide  modes  in  a conducting 


cylinder  of  radius  R^ . The  analysis  is  carried  out  for  a 


hollow 


e_  (3  C. 

electron  beam  (radius  R^)  propagating  parallel 


to  a 


uniform  axial  magnetic  field, xlt  is  assumed  that 

* U'vzN  \ 


v/y0«<L,  where  v is  Budker's  parameter  and  y^mc  is  the 


electron  energy  in  a frame  of  reference  moving  with  the  beam 


axial  velocity  B^ce^.  One  of  the  most  important  features 


of  the  analysis  is  that  the  instability  growth  rate  for 
magnetic  harmonic  numbers  s=2,3,...  is  comparable  to  the 
growth  rate  of  the  fundamental  (s=l)  mode,  particularly 
for  moderate  electron  energy  -(-P' 


■H — l /y“)  ^ ^>0. 41 . Moreover, 


0 v ' '0 

it  is  shown  that  the  instability  growth  rate  can  be  maximized  / 


by  appropriate  choice  of  the  geometric  parameter  R^/R^. 
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I.  INTRODUCTION 


One  of  the  most  basic  instabilities  that  characterizes  a hollow 

electron  beam^  with  slow  rotational  equilibrium  is  the  electron  cyclotron 

. „ . 2-10  _ J 4,5  , 

maser  instability.  In  recent  experiments,  the  cyclotron  maser 

instability  has  been  investigated  for  magnetic  harmonic  number 

s>2 , with  particular  emphasis  on  the  implications  for  intense  microwave 

generation.  Previous  theoretical  analyses  of  this  instability  have 

been  carried  out  for  a non-self-consistent  equilibrium  slab  configuration. 

Strictly  speaking,  a more  accurate  theoretical  model  of  microwave 

generation  by  the  cyclotron  maser  instability,  including  a determination 

of  the  optimum  value  of  the  beam  radius  R^,  requires  a linear  stability 

analysis  for  perturbations  about  a self-consistent,  cylindrical,  Vlasov 

equi 1 ibr ium. 

This  paper  develops  a self-consistent  theory  ^ of  the  cyclotron 

maser  instability  for  azimuthally  symmetric  perturbations  with 

magnetic  harmonic  number  s>l.  The  present  work  extends  the  previous 

self-consistent  theory^  of  the  cyclotron  maser  instability,  developed 

by  the  authors  for  s=l,  to  higher  values  of  magnetic  harmonic  number 

(s>2) . The  analysis  is  carried  out  within  the  framework  of  the  Vlasov- 

Maxwell  equations  for  an  infinitely  long  electron  beam  propagating 

parallel  to  a uniform  magnetic  field  B^e^  with  axial  velocity  V^e^. 

We  assume  that  the  beam  is  very  tenuous,  so  that  the  perturbed  field 

can  be  approximated  by  the  vacuum  waveguide  fields  where  v 

is  Budker's  parameter).  Equilibrium  and  stability  properties  are 

calculated  for  the  specific  choice  of  equilibrium  electron  distribution 

function  in  which  all  electrons  have  the  same  value  of  canonical 

2 

angular  momentum  (Pq)  and  the  same  value  of  energy  (yqIiic  ) in  a frame 


6,8 
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of  reference  moving  with  axial  velocity  fEq.  (5)].  Equilibrium 

properties  are  examined  in  Sec.  II.  An  important  feature  of  the  analysis 
is  that  the  equilibrium  distribution  function  in  Eq.  (5)  corresponds 
to  a hollow  density  profile  with  sharp  radial  boundaries  [Eq.  (13)]. 

The  formal  stability  analysis  for  azimuthal ly  symmetric  electro- 
magnetic perturbations  (;i/;)0=0)  is  carried  out  in  Sec.  III.  Making 
use  of  the  fact  tiiat  tiie  electron  trajectories  are  circular,  the  per- 
turbed distribution  function  is  calculated  for  arbitrary  magnetic 
iiarmonic  numbers  [Eqs.  (38)  and  (39)].  Equations  (38)  and  (39), 
when  combined  with  Eqs.  (16),  (22),  (23),  and  (24),  constitute  one 
of  the  main  results  of  this  paper  and  can  be  used  to  investigate 
stability  properties  for  a broad  range  of  system  parameters.  In 
this  regard,  we  emphasize  that  Eqs.  (38)  and  (39)  are  derived  with  no 

a priori  assumption  that  the  electron  motion  in  the  beam  frame  is 

2 2 

much  smaller  than  the  speed  of  light  in  vacuo,  i.e.,  that  6q  ;(1-1/yq)«1. 

In  Sec.  IV,  a detailed  analytic  investigation  of  the  cyclotron 

maser  instability  is  carried  out  for  the  TE  and  TM  wavegcide  modes, 

2 

assuming  6q<<1.  Introducing  the  normalized  Doppler  shifted  eigenf requency 
[Eq.  (49)] 

x =[y  (co-kV.  )— so)  ]/u)  , 

s b b c c 

the  dispersion  relation  can  be  expressed  as  [Eqs.  (48)  and  (51)] 

VQsnXs+e0sQsn/(2s+1)=0’  TE  mode’ 

Xl“^nxl+eoQlr/5=0’  ™ mode’ 

2 1/2 

where  y,=l/(l-B,  ) , w =eB»/Y„mc  is  the  electron  cyclotron  frequency 

b b c 0 0 

in  a frame  of  reference  moving  with  axial  velocity  S^ce^,  k is  the  axial 

L 


4 
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wavenumber,  gj  is  the  complex  eigenf requency , (Q^)  •*-s  t^le  coupling 
coefficient  for  the  TE  ( TM)  mode  [Eqs.  (50)  and  (52)],  and  s=l , 2,  3... 
denotes  the  magnetic  harmonic  number.  Evidently,  for  given  s which 


satisfies  a , <a_  for  the  TE  mode  and  a , <6„  for  the  TM  mode,  the 
si  On  si  On 

maximum  growth  rate  occurs  at  a value  of  R^/R^  given  by  [Eq.  (53)] 


Rn/R  = 
0 c 


a ,/a„  , TE  node  , 
si  On 

asl/60n’  ™ mode  ’ 


where  a . is  the  first  zero  of  J'(y)=0,  and  a_  and  are  the  nth 
si  On  On 

zeroes  of  J^(y)=0  and  ^q(y)=0>  respectively. 

A detailed  numerical  analysis  of  the  dispersion  relations  in  Eqs. 

(48)  and  (51)  is  presented  in  Sec.  V,  where  stability  properties  are 
investigated  for  a broad  range  of  system  parameters.  It  is  found 
that  the  growth  rate  of  perturbations  with  s>2  increases  rapidly 
when  the  value  of  (3^  is  increased.  We  therefore  conclude,  for  moderate 
or  high  values  of  (3  , that  magnetic  harmonic  perturbations  with  s>2 
are  also  important  unstable  modes  for  intense  microwave  generation  by 
the  cyclotron  maser  instability. 

Finally,  we  note  that  Eqs.  (53)  and  (30)  can  be  combined  to  determine 

the  conditions  for  maximum  growth  rate  and  hence  optimum  conditions 

for  intense  microwave  generation  by  the  cyclotron  maser  instability. 

In  particular,  the  microwave  frequency  io-sw^y^  produced  by  the  electron 

beam  can  be  tuned  by  the  matching  condition  R^ag^c/sto^  in  Eq.  (30). 

Moreover,  selecting  R./R  =a  , /a.  maximizes  the  growth  rate  of  the  TE 

U c si  On 

mode  perturbation  with  magnetic  harmonic  number  s (Sec.  IV  and  V). 
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I I . EQUILIBRIUM  PROPERTIES  AMI)  BASIC  ASSUMPTIONS 

As  illustrated  in  !•'  i >*, . I,  the  t‘<|u  i 1 i hr  ium  configuration  consists 
of  a slowly  rotating  hollow  electron  beam  that  propagates  parallel  to 
a uniform  applied  magnetic  field  with  mean  axial  velocity  3^c^z. 

The  mean  radius  of  the  electron  beam  is  denoted  by  and  a grounded 
cylindrical  conducting  wall  is  located  at  radius  r=R^.  The  applied 
magnetic  field  provides  radial  confinement  of  the  electrons,  and  the 
radial  thickness  of  the  electron  beam  is  denoted  by  2a.  As  shown  in 
Fig.  1,  we  introduce  a cylindrical  polar  coordinate  system  (r,0,z) 
with  the  z-axis  coinciding  with  the  axis  of  symmetry;  r is  the  radial 
distance  from  the  z-axis,  and  9 is  the  polar  angle  in  a plane  perpendicular 
to  the  z-axis. 

The  following  are  the  main  assumptions  pertaining  to  the  present 
analysis : 

(a)  Equilibrium  properties  are  independent  of  z (3/3z=0)  and 
azimuthally  symmetric  (3/90=0)  about  the  z-axis. 

(b)  The  mean  canonical  angular  momentum  of  the  electrons  is 
negative,  which  corresponds  to  a slow  rotational  equilibrium. ^ ^ 

(c.)  It  is  further  assumed  that 


v/y0<«l  (1) 

2 2 

where  v=N^e  /me  is  Budker’s  parameter, 

N =2tt  f c dr  r n°(r)  (2) 

e J0  e 

is  the  number  of  electrons  per  unit  axial  length,  n^(r)  is  the 

e 

equilibrium  electron  density,  c is  the  speed  of  light  in  vacuo, 

-e  and  m are  the  charge  and  rest  mass,  respectively,  of  an  electron. 


and  y^mc.  is  the  electron  energy  in  a frame  of  reference  moving  with 
the  mean  axial  velocity  of  the  electron  beam.  The  inequality 

in  Eq.  (1)  indicates  that  the  beam  is  very  tenuous,  so  that  the 
perturbed  fields  can  be  approximated  by  the  vacuum  waveguide  fields."^ 
Consistent  with  the  low-density  assumption  in  Eq.  (1),  we  also  neglect 


the  influence  of  the  small  equilibrium  self-electric  and  self-magnetic 
fields  that  are  produced  by  the  lack  of  equilibrium  charge  and  current 
neutralization. 

For  azimuthally  symmetric  equilibria  with  3/30=O=3/3z,  there  are 
three  single-particle  constants  of  the  motion.  These  are  the  total 
energy  H, 

H=Ymc  =(mc+cj))  , (3) 


the  canonical  angular  momentum  P , 

D 


Pe==r  t Pq  Ce/2c)  rB0 1 


and  the  axial  canonical  momentum  P =p  . In  Eqs.  (3)  and  (4),  lower 

z z 

case  £ denotes  mechanical  momentum  and  the  equilibrium  self-fields 
have  been  neglected  in  comparison  with  the  external  magnetic  field 
[see  Eq.  (1)],  Any  distribution  function  that  is  a function 
only  of  the  single-particle  constants  of  the  motion  satisfies  the 
steady-state  Vlasov  equation  (3/3t=0).  For  present  purposes, 
we  assume  an  equilibrium  distribution  function  of  the  form,^’'*'^ 

f°(».p9.pz)  - Sft»«w  , (5) 

2“Vo 

2 2 

where  nQ=const.  is  the  electron  density  at  r=RQ.  Pg=-(e/2c) (R^-a  )Bg=const 
is  the  canonical  angular  momentum  of  the  electrons, 


u=H"ebcVYomc  /Yb 
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is  an  effective  energy  variable,  B =const.  is  defined  by  3 =V,  /c= 

n b b 

2 l <2 

(y, -1)  /y,  , and  y„  and  y,  are  constants.  It  is  straightforward 

b b u 1) 

to  show  that  the  axial  velocity  profile  associated  with  Eq.  (5)  is 


uniform  over  the  beam  cross  section. 

Several  pertinent  equilibrium  properties  can  be  deduced  for 

the  class  of  electron  beam  equilibria  described  by  Eq . (5).  For 

this  purpose,  it  is  useful  to  transform  the  energy  variable  U 

defined  in  Eq.  (6)  from  momentum  variables  (p  ,pg,p  ) appropriate 

to  the  laboratory  frame  to  momentum  variables  (p',p',p')  appropriate 

r 0 z 

to  a frfmie  of  reference  moving  with  velocity  3,  ce  , where  e is  a 

J b 'vz  'vz 

unit  vector  in  the  z-direction.  The  relevant  transformation^’^ 
is  given  by 


Pr=Pr’  p0=p<r  Pz=>b(pz+Y’m6bc)  * Y=VY ’+6bPz/mc)  , (7) 


U = ~ (y'-yn)mc2  , (8) 

Yb  ° 

2221/2  2221/2 

where  y=(l+j>  /me)  and  y'  = (l+p'  /me)  . After  some  straight- 
forward algebra,  we  find 

n ? n Rn(y'+6  p'/mc)  ? ? a 

fe(H’P0’Pz)d  P = ~iimY0Y’ 6(Y'mc  “Y0mC  )6(P(TP0)d  p'  ’ 

(9) 

3 3 

where  d p=dpfdpgdpz  and  d p ' =dpMp^dp^ . It  is  evident  from  Eq.  (9) 

2 

that  YqItic  can  be  identified  with  the  total  electron  energy  in 

a frame  of  reference  moving  with  axial  velocity  8^c. 

2 

The  energy  variable  y'mc  in  Eq . (9)  can  be  expressed  as 


, 2 . 2 2X  2-2  2 ,2,1/2 

y me  =c(m  c +m  +p|  ) 


where  use  has  been  made  of  P^=Pq,  g is  defined  by 


8 


£ 

g(r)  = -°-  [p+(p2+a2)/2RQ]  , 


(ID 


wc=eB()/mc  is  the  nonrelat iv istic  electron  cyclotron  frequency, 
and  p|2=p^2+p^2.  In  Eq.  (11)  p is  defined  by 


p=r-R0  . 


(12) 


Substituting  Eq.  (11.)  into  Eq.  (9),  and  representing 


d3p'  = 


2u 

dPg | da  I dpjp|  , 


'0  J0 

it  is  straightforward  to  show  that  the  electron  density  can  be 
expressed  as 

f,3  ,0.  R- 


0,  , 

ne(r)= 


where 


d-pr(H,Pe,Pz)  = nQ  «(a2-p2)  , 


, 2 .,1/2  r 
a=(y  -1)  c/w 
U c 


(13) 


(14) 


is  the  half-thickness  of  the  beam,  and  «(x)  is  the  Heaviside  step 
function  defined  by 

( 0 , x<0  , 

• = 

v 1 , x>0  . 

The  self-consistent  electron  density  profile  is  illustrated  in  Fig.  2. 
Evidently,  the  electron  beam  equilibrium  described  by  Eq.  (5)  has 
sharp  radial  boundaries.  Additional  equilibrium  properties  associated 
with  the  distribution  function  in  Eq.  (5)  are  discussed  in  Refs.  1 


and  10. 
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1 1 1 . LINEARIZED  VLASOV-NAXWELL  EQUATIONS  FOR  TENUOUS  BEAM 


Id  this  section,  we  make  use  of  the  linearized  Vlasov-Maxwei 1 
equations  to  investigate  stability  properties  for  azimuthally  symmetric 
perturbations  (3/30=0)  about  a tenuous,  hollow-beam  equilibrium 
described  by  Eq.  (5).  We  adopt  a normal-mode  approach  in  which  all 
perturbations  are  assumed  to  vary  with  time  and  z according  to 


6^(x,t)=i(j(r)exp{i(kz-wt) } , 


where  lmu>>0.  Here  w is  the  complex  eigenf requency  and  k is  the  axial 
wavenumber.  The  Maxwell  equations  for  the  perturbed  electric  and  magnetic 
field  amplitudes  can  be  expressed  as 


VxB(x)  = — J(x)-i  - E(x) 
% % % c 'v  ^ c % 


where 


J(x)=-e  d P^fe(x,£) 


is  the  perturbed  current  density.  In  Eq . (16), 


diexp{-icoi  }<  £(xf ) + 


8 0 

3^’  e 


is  the  perturbed  distribution  function,  T=t'-t,  and  the  particle 
trajectories  x'(t')  and  j^'(t')  satisfy  d^'/dt'=v'  and  dj^'  /dt ' =-ev' xB^e^/c , 
with  "initial"  conditions  x'(t'=t)=x  and  v'(t'=t)=v. 

In  the  tenuous  beam  limit  consistent  with  Eq.  (1),  the  perturbed 
fields  can  be  approximated  by  the  vacuum  waveguide  fields.  ^ 

In  this  context,  the  present  stability  analysis  utilizes  the  vacuum 
transverse  electric  (TE)  and  transverse  magnetic  (TM)  waveguide  modes 


10 


as  a convenient  basis  to  represent  a general  electromagnetic  field 
perturbation  within  a cylindrical  waveguide.  Making  use  of  Eq.  (15) 
and  neglecting  the  perturbed  current  density,  the  vacuum  waveguide  modes 
can  be  expressed  as 


Bz(r)=Jo(aonr/Rc)=-i(c/wrH3(rEe)/8r]=i[3(rBr)/8r]/kr  , (18) 


for  the  TE  mode  and 


Ez(r)=Jo^Onr/,Rc^  = i^c/,<J)r^3^rBe^8r^=i^3^rEr^9r^kr  ’ (19) 

for  the  TM  mode.  In  Eqs.  (18)  and  (19),  J^(x)  is  the  Bessel  function 

of  first  kind  of  order  £,  and  a.  and  8_  are  the  nth  roots  of  J, (a„  )=0 

Un  On  1 On 

and  respectively.  Moreover,  without  loss  of  generality, 

the  normalization  amplitudes  for  Bz(r=0)  and  Ez(r=0)  have  been  set 
equal  to  unity  in  Eqs.  (18)  and  (19).  After  a simple  algebraic 
manipulation  of  Eqs.  (15),  (18),  and  (19),  it  is  straightforward 


to  show  ti  at 


^ ’ IF)  Jl(aOnr/Rc)="47T(aOn/Rcc)J0(r)  ’ (20) 

C 

for  the  TE  mode  and 

~ k2  " ^r1)  J0(B0nr/Rc)=47rik[pe(r)-(u./kc2)Jz(r)]  (21) 

\ c R^  / C 

c 

for  the  TM  mode.  In  Eq . (21),  the  perturbed  charge  densitv  p (r)  is 

e 


defined  by 


Pe(r)=-e[d3pfe(r,£) 


Multiplying  Eqs.  (20)  and  (21)  by  rJ, (a  r/R  ) and  rJn(3_,  r/R  ), 

l un  c u Un  c 

respectively,  and  integrating  from  r=0  to  r=R  , we  obtain 

c 
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J .2  U0n' 

“2  “ k " ~ 
C R ) 

c 


for  the  TE  mode  and 

f uj2  2 p0n' 

\ c R / 

c 


8iia_  /R  c 
On  c 


[R(,.2(a0n)]  J0 


' dr  r Jl(aOnr/Rc)J0(r)  (23) 


8ti  ik 


R 


, I 1 dr  r J (8„  r/R  ) p (r)  J 

[R  j.(e  )]2  Jo  0 0n  c L e kc2 

c 1 On  (24) 


for  the  TM  mode.  For  present  purposes  it  is  also  assumed  that 


I ft  I = I w-kV.  -sw  /y,  I <<u>  /y, 
s b c b c b 


(25) 


where  ft  “w-kV^-sw  /y^  is  the  Doppler  shifted  frequency,  wc=eB^/yQmc 
is  the  electron  cyclotron  frequency  in  a frame  of  reference  moving 
with  axial  velocity  V^,  and  s=l,2,3...  denotes  the  magnetic  harmonic 
number. 

To  lowest  order,  the  eigenf requency  w and  axial  wavenumber  k 
are  obtained  from  the  simultaneous  solution  of  the  vacuum  waveguide 
mode  dispersion  relation, 

2 


w , 2 

~2  “ k 
c 


2 2 

| a0n^Rc  » TE  mode  > 

02  /R2  , TM  mode  , 
On  c 


(26) 


and  the  condition  for  cyclotron  resonance 


w^kV, +sw  /y,  . 

b c b 


(27) 


Moreover,  to  maximize  the  growth  rate  and  efficiency  of  microwave 
generating  it  is  required  that  the  group  velocity  of  the  vacuum 
waveguide  mode  in  Eq.  (26)  be  approximately  equal  to  the  beam 
velocity,  i.  e. , 


V = % = - V, 

g dk  u b 


(28) 
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Solving  Eqs.  (27)  and  (28)  for  the  characteristic  frequency  and  axial  wave 
number  (w,  k)  = (w^  , k^) , we  find  (Fig.  3) 


VsVb 


ko=sVb&b/c 


(29) 


2 -1/2 

where  y^=(l-S^)  . For  maximum  growth,  it  is  also  required  that 

(iOq^q)  solve  Eq.  (26)  in  leading  order.  Therefore,  for  maximum 
growth,  we  find  that  should  satisfy 

!an  c/sto  , TE  mode  , 

0n  C (30) 

c/sto  , TM  mode  , 

On  c 


for  intense  microwave  generation  at  frequency  u^sto^Y^-  Because  of  the 
discrete  nature  of  a..  , , and  s,  we  also  conclude  that  the  stability 

analysis  can  be  carried  out  separately  for  the  TE  and  TM  modes. 

The  perturbed  distribution  function  is  calculated  for  the  case  of 
a self-consistent  Vlasov  equilibrium  in  which  all  electrons  have  the 
same  canonical  angular  momentum  and  the  same  total  energy  in  a frame  of 
reference  moving  with  axial  velocity  V^z  [Eqs.  (5)  and  (9)].  To 
simplify  the  right-hand  side  of  Eq.  (17),  we  make  use  of  Eqs.  (4) 
and  (6),  and  the  identities  3U/8p=v-V.  e and  3P./3p^=re„,  where  e„ 
is  a unit  vector  in  the  ©-direction.  The  TE  mode  portion  of  the 
perturbed  distribution  function  can  then  be  expressed  as"^ 


(31) 


which  is  required  to  calculate  the  perturbed  azimuthal  current 
density  Jg(r)  in  Eq-  (23).  Similarly,  the  TM  portion  of  the  perturbed 
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distribution  function  is  given  by 

iU 


10 


3f°  r 0 

f“(r,£)=e  • drexp{i[(kp  /ym)-o)]T}{[Er(r,)-6bB0(r')]vJ. 

J — a> 


(32) 


+ Ez(r')  t (pz / Yni) ~Vb ] } , 


which  is  required  to  calculate  the  perturbed  charge  and  axial  current 
densities  in  Eq.  (24).  In  obtaining  Eqs.  (31)  and  (32),  use  has  been 
made  of  the  axial  orbit 


z'=z+(p7/-ym)  ( t ' -t)  . 


The  transverse  (radial  and  azimuthal)  motion  of  a typical  electron 
is  illustrated  in  Fig.  4.  (The  dotted  circle  is  the  electron  orbit 
in  a plane  perpendicular  to  the  z-axis.)  The  radial  distance  of  the 
electron  from  the  z-axis  at  times  t'=t  and  t'=t'  are  denoted  by  r and 
r',  respectively.  The  point  C is  the  gyrocenter  of  the  electron 
trajectory.  The  angular  coordinates  (j>  and  <J) ' are  the  perpendicular 
velocity-space  polar  angles  at  times  t'=t  and  t'=t',  and  are  related  by 


>’=<t>+(uc/y)  (t'-t) 


The  transverse  velocities  at  times  t'=t  and  t'=t'  are  denoted, 

respectively,  by  v^,  and  v^, , and  the  corresponding  speeds  are  defined  by 

. 2 2.1/2  , , , ,2,  ,2,1/2 

VT=(W  and  v;=(v;  +v^  ) . 

To  simplify  the  present  analysis,  we  also  assume  that 


v/V<(B0wcRc/c>‘ 


(33) 


2 1/2 

where  6 =(1-1/Yq)  . Equation  (33)  is  easily  satisfied  in  parameter 

2-5 

regimes  of  experimental  interest.  Within  the  context  of  Eq.  (33), 
it  is  valid  for  s<3  to  neglect  the  terms  proportional  to  9f^/3P  in 
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Eq.  (31),  since  the  corrections  associated  with  these  terms  are  of  order 
(v/Yg) (c/£qW  R ) “ (<<1)  or  smaller.  Substituting  Eq.  (18)  into 
Eq.  (31),  we  find  that  the  perturbed  TE  mode  distribution  function 
can  be  approximated  by 


ieR  p„  3f  rO 

r (r ,p)  = (co-kV.)  — dtexp{i[(kp  /ym)-w]T} 

e ^ c b ym  311  J z 

On  J -°° 


xsin(<{i ' -0 ' ) (aQnr ' /R£)  , 


221/9  f y 

where  p =ymv1j,=  (p*’+Pg)  and  use  has  been  made  of  Vg=v^,sin((ji ' -6  ). 

Making  use  of  Eq.  (28)  and  substituting  Eq.  (19)  into  Eq.  (32),  we 


obtain 


3f ° fO 

f”(r,^)=e[(pz/ym)-Vb]  — — j diexp{i [ (kp2/Ym)-u)]T}J0(60nr' /Rc) 


for  the  perturbed  TM  mode  distribution  function. 

The  Bessel  function  summation  theorem  for  the  triangle  OAB 
in  Fig.  4 can  be  expressed  as^ 

exp{U(0'-o)}J£(aOnr'/Rc)=  J J l+g  ' (a0nr/Rc)^ Js ' (a  OnX/V 

S =-oo 


xexp{is'  (n+9-<p-ucT /2y)  } 


Similarly,  for  the  triangle  ABC,  we  can  represent 
exp  { i (tt-u^t /y)  } , (a^A /Rc) 

(37 

= j J J exp{ism  x/y } , 

L s+s  Vmio  R / slmii)  R / c 

s=-<»  c c c c 

where  s=l,2,3...  denotes  the  magnetic  harmonic  number.  Making  use  of 
Eqs . (28),  (36),  and  (37),  and  carrying  out  the  time  integration 
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in  Eq.  (34),  it  is  straightforward 
component  s,  that  the  perturbed  TE 
imated  by 


eR  (w-kV,  ) p 3f° 
c b 1 e 

a„  c 2y,  m 3U 

On  b 


to  show,  for  a given  harmonic 

mode  distribution  function  can  be  approx- 


Wt\_  3 /Vt 

mo)  R / s+1  V mu)  R 

c c c c 

Y1  (u-kVb)-s£i  /Yb 


) 


(38) 


^,(i)S'Js’(aOnr/Rc)Js+s'  (^|1r)eXp{is,(e^)}  * 
s c c 

where  use  has  been  made  of  the  Lorentz  transformation  in  Eq.  (7). 

In  obtaining  Eq.  (38),  we  have  neglected  the  mode  coupling  between 
different  values  of  s,  which  is  consistent  with  Eq.  (25).  In  a 
similar  manner,  the  perturbed  TM  mode  distribution  function  for  harmonic 
component  s can  be  approximated  by 


Equations  (38)  and  (39),  when  combined  with  Eqs.  (16),  (22), 

(23),  and  (24),  constitute  one  of  the  main  results  of  this  paper 

and  can  be  used  to  investigate  stability  properties  for  a broad  range 

of  system  parameters.  In  this  regard,  we  emphasize  that  Eqs.  (38) 

2 2 

and  (39)  have  been  derived  with  no  a priori  assumption  that  (1-1/Yq) <<1 

2 

or  that  a<<Rg.  However,  in  the  limiting  case  where  6q<<1,  Eqs.  (38) 
and  (39)  can  be  simplified  considerably  (Sec.  IV). 
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IV.  CYCLOTRON  MASER  INSTABILITY  KOR  Eq<<1 

In  this  section,  simplified  expressions  are  obtained  for  the 

perturbed  distribution  function  in  Eqs.  (38)  and  (39),  and  the 

results  are  used  to  derive  the  dispersion  relation  for  several  different 

values  of  magnetic  harmonic  numbers.  The  present  analysis  assumes 

that  the  electron  motion  in  a frame  of  reference  moving  with  axial 

velocity  B c is  much  smaller  than  speed  of  light  in  vacuo,  i.e.,  that 
b 

Bq=1-1/Yq«1  . (AO) 

Equation  (40)  can  be  used  to  truncate  the  summations  ov< r s'  in 

2 

Eqs.  (38)  and  (39),  keeping  only  leading  terms  of  order  Bq- 
E M 

Defining  Q and  Q to  be  the  coupling  coefficients  between  the 
sn  sn  ° 

vacuum  TE  and  TM  waveguide  modes  and  the  cyclotron  resonance  mode 

w=kV,+su  /y  , it  is  found  that  |qM  /Q^‘  I <6^/4  [see  Eqs.  (50)  and  (52)]. 
b c b sn  sn  “v  U 

2 

In  this  context  (Bq<<1),  the  instability  growth  rate  is  largest 
for  the  TE  mode. 

The  perturbed  distribution  function  in  Eqs.  (38)  and  (39)  can 

be  further  simplified  by  making  use  of  the  symmetry  properties  of 

the  equilibrium  distribution  function  f^(U,P„).  Since  the  variable 

e 0 

U is  an  even  function  of  p^  [Eq.  (8)],  it  follows  from  Eqs.  (16) 
and  (22)  that  any  term  in  Eqs.  (38)  and  (39)  that  is  an  odd  function 
of  p^  will  give  zero  when  the  integration  over  p^  is  carried  out. 
Therefore,  when  evaluating  fes(r>£)»  we  simply  omit  terms  proportional 
to  odd  functions  of  p^.  To  evaluate  the  momentum  integral  in  Eqs.  (16) 
and  (22),  use  is  made  of  Eq . (10).  After  some  straightforward  algebra. 


we  f ind 
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,3  1 

d p 


3f°/3U 

e 


noRo 


y y * (o>-kV,  ) -sdj  /y. 

b cb 


+ (l  - -— ) <5 (p+a) 


W 


an 


(l  + -~)  6 (p-a) 


o P 
3 rr 


d p 


0 

3f°/9U 

e 


VW“kV  ,2  2. 

— ^ © (a  -p  ) 

afi 

S 


fa3  r 


,2 


y' (w-kV,  )-su  /y,  'b 


3 f ° / 9U 
e 


d P y'  y'  (o)-kVb)-swc/Yb 


d3  Pr 

d P T 


n0R0m 

V 

9f°/3U 

e 


g2(a2-g2)  (w-kVb)'| 


2 2 
2a  n 

s 


2 2 

©(a  -p  ; , 


Y y' (w-kVb)-swc/Yb  Yb 
3 2 

3n0R0m  “c  , 2 2/ 

27^ — (a  "8  } 

Y0 


= — d P 


,3  , <!»2 


9f°/3U 

e 


y'  (w-kV.  )-sai  /y, 
b c b 


e2(w-kVb)(l-g2/a2) 


(41) 


4n 


2 2 

©(a  -p  j , 


where  the  Doppler  shifted  frequency  ng  is  defined  in  Eq.  (25),  and 
g(r)  is  defined  in  Eq.  (11). 

We  first  investigate  stability  properties  for  the  TE  mode  pertur- 
bation with  fundamental  magnetic  harmonic  number  (s=l) . From  Eq.  (38), 
the  perturbed  distribution  function  is  approximated  by 


'fel 


eR  w-kV,_ 
c h 


f".  (r ,p)  = - — — ~ 11  3,  (u  r/R  ) — 

^ a c 2y,m  1 On  " - 


P8Sf"/30 


On 


c-  Y'(a)-kVb)-<WYb  ’ 


(42) 


where  pf)=p^,sin(4)-0)=mwcg.  The  perturbed  azimuthal  current  density 
is  evaluated  by  substituting  Eq . (42)  into  Eq . (16)  and  making  use 
of  Eq.  (41).  After  some  algebra,  we  obtain 


E . . n0e  R0Rc(a)~kVb) 
J91U;  " 2YnYuman  cr 


J.(«n  r/R Jg2  ]-4-  (l  + 6(p-a) 

0'b““0n“  1 °n  C f3“l  ^ V 


(43) 


+ (l  - — ) 6 (p+a) 


B0(u"kVb)  , 2 2. 

— j~2 ©(a  ) 

a fi 


_ 
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The  radial  integration  on  the  right-hand  side  of  Eq.  (23)  can  be 

carried  out  by  Taylor  expanding  the  integrand  about  r=RQ,  and  keeping 

2 2 

leading  terms  of  order  p =(r-Rn)  . The  Taylor  expansion  of  the  Bessel 
functions  provides  a good  approximation  when  the  functions  are  slowly 
varying  over  the  minor  cross  section  of  the  beam.  The  number  of 
electrons  per  unit  axial  length  of  the  beam  (N^)  can  be  determined 
by  substituting  Eq.  (13)  into  Eq.  (2),  which  gives  N^^Anan^Rg. 
Eliminating  n^  in  favor  of  and  carrying  out  some  straightforward 
algebraic  manipulation,  we  obtain  the  s=l  TE  mode  dispersion  relation 


? z 
u ,2  a0n  2v 

fV'WVi 

2 

x>\ 

> 

3 

c2  R2  'Vb 

'RcJ2(a0n)  J 

Ico-kV, -aj  /y, 

L b c b 

c 


6o(“-kV2 


(44) 


3(w-kV  -uj  /y,)2-* 
b c b 

2 2 

where  use  has  been  made  of  v=Nee^/mc  . The  dispersion  relation  in 
Eq.  (44)  is  identical  to  the  result  obtained  previously  by  Uhm  et  al. 
for  s=l . 


In  a similar  manner,  we  have  derived  the  dispersion  relation 

for  several  values  of  magnetic  harmonic  number  s.  The  following 

2_  2 

equations  summarize  the  main  analytic  results  obtained  for  Bq=(1-1/yq) <<1 
(a)  The  TM  mode  dispersion  relation  for  s=l  is  given  by 


u>2  .2  B0n 

t - k - ir 

c R 

c 


2vB 


0 


, 3 

3Vb 


Jl<60nVRc> 


Vl(V 


“-kV“c/,b 


8Qw(w-kVb) 


(45) 


5(a,-kVb-(oc/yb)  J 

2 

where  use  has  been  made  of  the  recursion  formula  J^-d( J^J^) /dx- 


Vi/X-Ji- 
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(b)  The  TE  mode  dispersion  relation  for  s=2  is  given  by 


SL‘  - k2  - - 

c2  R2 

c 


On 


8'"’o 

'’o'1!, 


Rc82(a0n) 


6o<--kV2 

5(w-kVb-2wc/Yb)  ; 


where  use  has  been  made  of 


w-kV, 


w-kV, -2u>  /y, 
b c b 


J2-d(J1J2)/dx-J1J2/x=J2 


(46) 


(c)  The  TE  mode  dispersion  relation  for  s=3  is  given  by 


- k‘ 


T)n 
2 " 


81v3q_ 

2°Vb 


J3(ra0nR0/Rc) 


R J„ (an  ) 
c 2 On 


to-kV, 


ui-kV  -3w  jfy, 

b c b 


Bo(“-kV2 


7(w-kVb-3wc/Yb)' 
where  use  has  been  made  of 


(47) 


J2=J2-2d  (J  J ) /dx-2J1  J2/x+d2(J1  J3)  /dx2+3  [d(J1J3)  /dx]  /x  . 

Substituting  Eq.  (30)  into  Eqs.  (44),  (46),  and  (47),  and  making 

use  of  u)-kV, +sw  / y,  , we  obtain  the  approximate  TE  mode  dispersion 
b c b 

relation 


x3-QE  x +e„sQE  / (2s+l)  = 0 
s sn  s 0 sn 


(48) 


for  s=l,2,3.  In  Eq.  (48),  the  normalized  Doppler-shifted  eigenf requency 
xg  is  defined  by 


V[Yb(u)-kV-SUc]/“c  ’ 


(49) 


and  the  coupling  coefficients  between  the  TE  vacuum  waveguide  mode 
and  the  electron  cyclotron  resonance  mode  are  defined  by 
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Qln(R0/Rc)  V^Jl(ci0nR0/Rc)/J2(a0n)  1 /Y0Yba0n  ’ 


4(R0/Rc)  = 16^0[J2(a0nR0/Rc)/J2(a0n)]2/3W0n  * (50) 


Q3'n(R0/Rc)  9vO60)  tJ3(°lOnRO^Rc^J2^a0n^  ^0Y0Yb“0n 


Similarly,  for  the  s=l  TM  mode,  we  obtain 


1 M 2 M 

XrQlnXl+Wn/5=°  • 


(51) 


where  the  coupling  coefficient  is  defined  by 


QL(R0/Rc)-''B0[Jl(60nE0/V/Jl<60n>>2/3Vb60„ 


(52) 


As  mentioned  at  the  beginning  of  this  section,  it  is  evident 

from  Eqs.  (50)  and  (52)  that  the  TM  mode  coupling  coefficient 

M E 

is  much  less  than  the  TE  mode  coupling  coefficient  Q^. 

Therefore,  we  conclude  that  the  TE  mode  is  the  dominant  unstable  pertur- 
2 

bation  for  Bq<<1.  However,  when  the  transverse  electron  motion  is 

relativistic  (Sq^I),  the  TM  mode  perturbation  is  equally  important. 

A careful  examination  of  Eqs.  (50)  and  (52)  shows  that  the  coupling 

coefficient  between  the  vacuum  waveguide  mode  and  the  electron 

cyclotron  resonance  mode  (w=kV^+  scWy^)  is  a maximum  whenever 

J'(a„  R„/R  )=0  for  the  TE  mode,  and  J ’ ( B „ R„/R  )=0  for  the  TM  mode, 
s On  0 c s On  0 c 

Here  the  prime  (')  denotes  dJs(x)/dx.  In  this  context,  we  find  that 

the  maximum  growth  rate  for  magnetic  harmonic  number  s occurs  for  a 

value  of  R_/R  given  by 
0 c 


Rn/R 
0 c 


sl^On  ’ 

TE 

mode 

sl/60n  ’ 

TM 

mode 

(53) 


where  is  the  first  root  of  J^(as3)=0.  Equation  (53)  is  valid  only 
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r v 


when  a , <a_  for  the  TE  mode,  and  a ,<P„  for  the  TM  mode.  For  a , >a„ 
si  On  si  On  si  On 

(TF.)  or  a (TM) , the  maximum  growth  rate  occurs  for  R„/R  =1. 

s 1 On  U r 

It  should  also  be  noted  from  Eq.  (50)  that  perturbations  with  higher 
magnetic  harmonic  number  become  dominant  when  the  transverse  electron 
speed  approaches  c (Sq^I).  In  the  following  section,  we  make  use  of 
Eqs.  (48)-(52)  to  investigate  detailed  stability  properties  for  a 
broad  parameter  range  of  experimental  interest. 
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V . NUMERICAL  ANALYSIS  OF  STABILITY  PROPERTIES 

The  growth  rate  w.=lmw  and  real  oscillation  frequency  wr=Rem 
have  been  calculated  numerically  from  Eqs.  (48)  and  (51)  for  a broad 
ange  of  system  parameters,  R^/R^,  n,  s,  BQ,  and  8^.  Since  the 
electron  beam  is  hollow  (R^>a)  and  is  located  inside  a conducting 
waveguide  of  radius  (R^R^-a),  the  allowable  range  of  Rq/Rc  is 
restricted  to 

(sg0/ci0n)<(R0/Rc)<1"sS0/a0n’  TE  mode  ’ 

(54) 

(sW<<E0/Rc)<1-s6O/6On'  ™“de  ’ 

where  use  has  been  made  of  Eq.  (30).  Therefore,  in  Figs.  5-7,  the 

plots  are  presented  only  for  values  of  R^/R  satisfying  Eq.  (54). 

Shown  in  Fig.  5 are.  plots  of  (a)  the  normalized  growth  rate 

and  (b)  the  normalized  Doppler-shifted  real  frequency 

Re(Y.fi,/u  )=[Yi_(w  -kV,  )-w  ]/w  versus  R„/R  obtained  from  Eq.  (48) 
blc  b r b c c 0c 

(TE  mode)  for  s=l,  v=0.001,  6 =0. 4,  6,  =0.286,  and  severa  values  of  n. 

U b 

Several  features  are  noteworthy  from  Fig.  5.  First,  the  maximum 

growth  rate  decreases  as  the  radial  harmonic  number  n is  increased. 

As  evident  from  Eqs.  (50)  and  (52),  this  feature  represents  a general 

tendency  for  both  TE  and  TM  perturbations,  for  all  magnetic  harmonic 

numbers  s.  Second,  the  number  of  zeroes  of  ok  (where  the  growth  rate 

vanishes)  increases  as  n is  increased.  Moreover,  the  maximum  value  of 

the  growth  rate  is  also  a decreasing  function  of  R^/R( . Third,  the  plot 

of  the  Doppler-shifted  real  frequency  versus  R^/R  has  the  same  general 

form  as  the  plot  of  the  growth  rate  versus  R„/R  . However,  the 

U c 

growth  rate  is  about  1.5  times  as  large  as  the  Doppler-shifted 

real  frequency.  This  feature  is  a result  of  the  fact  that  the  contribu- 


tion  of  the  second  term  in  Eq.  (48)  is  negligible  in  comparison  with 

the  third  term,  for  8^=0. 4.  Fourth,  the  maximum  growth  rate  occurs  at 

a value  of  R^/R^  corresponding  to  the  maximum  value  of  the  perturbed 

azimuthal  electric  field  [Eqs.  (18)  and  (44)). 

Shown  in  Fig.  6 are  plots  of  the  normalized  growth  rate  y^w^/o^ 

versus  R„/R  obtained  from  Eq.  (51)  (TM  mode)  for  parameters  identical 
0 c 

to  Fig.  5.  Comparing  Fig.  6 with  Fig.  5(a),  it  is  evident  that  the  TE 

mode  is  more  unstable  than  the  TM  mode.  For  example,  for  n=l,  the 

maximum  TE  mode  growth  rate  is  uk -0.0143  u^/y^  whereas  the  maximum 

TM  mode  growth  rate  is  u. -0.0055  oj  /y.  . In  this  context,  we  conclude 
° l c b 

that  the  TE  mode  is  a more  effective  means  for  exciting  microwave 

radiation,  at  least  for  modest  values  of  8^.  Since  the  plot  of  the 

Doppler-shifted  real  frequency  has  a similar  form  to  the  growth  rate 

curve,  it  is  not  shown  in  Fig.  6 (see  Fig.  5). 

We  now  present  an  example  that  illustrates  cyclotron  maser 

stability  properties  for  higher  magnetic  harmonic  numbers  (s>2) . 

Shown  in  Fig.  7 are  plots  of  the  normalized  growth  rate  y^w^/k^ 

versus  Rq/R  obtained  from  Eq . (48)  (TE  mode)  for  (a)  s=2  and 

(b)  s=3,  and  parameters  otherwise  identical  to  Fig.  5.  As  evident 

from  Fig.  7(a),  the  maximum  growth  rate  for  s=2  occurs  for  Rq/Rc=0.78 

when  n=l,  for  Rn/R  =0.44  when  n=2,  and  for  R„/R  =0.3  when  n=3.  These 
0 c 0 c 

values  of  R./R  correspond  to  R„/R  =a„-. /a_  [see  Eq.  (53)].  In  Fig. 

0 c 0 c 21  On 

7(b)  (s=3) , the  maximum  growth  rate  for  n=2  and  n=3  also  occur  at 
R„/R  =0.6  and  R~/R  =0.41,  respectively.  However,  for  s=3  and  n=l, 
the  growth  rate  assumes  a maximum  value  for  Rq/Rc=0.69,  since  a32>aop 
[see  Fig.  7(b)).  Evidently,  the  value  of  R^/R^a^/ oi^  plays  a very 
important  role  in  determining  optimum  system  parameters  for  intense 


microwave  generation  by  the  cyclotron  maser  instability. 
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Of  considerable  experimental  interest  is  the  stability  behavior 
for  specified  n and  several  values  of  s.  Typical  results  are  shown 
in  Fig.  8 where  (a)  the  normalized  maximum  growth  rate  ( y^uj ^ /a)c ) m 
and  (b)  the  normalized  Doppler-shifted  real  frequency  Re(-y^ft/wc) 
are  plotted  versus  8q  for  n=2  and  v=0.001.  In  Fig.  8,  we  assume 
Bg/1.4,  and  the  solid  and  dashed  curves  represent  the  TE  and  TM 
modes,  respectively.  The  range  of  8q  is  limited  to  0.1<Bq<0.5 

2 2 

since  this  stability  analysis  is  valid  only  when  (v/Yq) (c/R^oj^  <<8q<<1 

[see  Eqs.  (33)  and  (40)].  In  Fig.  8(b),  the  real  frequency  for  the 

s=l  TE  mode  is  related  only  for  the  range  of  Bg  corresponding  to 

instability  (uk>0) . For  the  TE  mode  perturbation,  maximum  growth 

occurs  at  R^/R  =0.26  for  s=l,  at  R„/R  =0.44  for  s=2,  and  at  R^/R  =0.6 
0 c 0 c 0 c 

for  s=3,  whereas  the  TM  mode  perturbation  with  s=l  has  a maximum 

growth  rate  at  Rg/Rc=0.33  [see  also  Figs.  5-7 J . Note  that  the 

s=l  TE  mode  perturbation  is  stabilized  by  decreasing  8q  to  Bg=0. 18, 

which  corresponds  to  Bg^^Q^/^) (see  Ref.  ]0).  As  shown  in  Fig.  8(a),  the 

s=l  TE  mode  perturbation  is  the  most  unstable  mode  for  small  values 

of  8q  (Bq<0.3,  say).  However,  the  growth  rate  of  perturbations  with 

s>2  increases  rapidly  when  the  value  of  Bg  is  increased.  As  an  example, 

for  the  (s,n)=(3,2)  TE  mode,  the  maximum  growth  rate  increases  by  a 

factor  of  six  when  Bq  is  increased  from  0.2  to  0.5.  We  therefore 

conclude  that  for  moderate  or  high  values  of  Bg,  perturbations  with 

higher  magnetic  harmonic  numbers  (s>2)  are  also  important  unstable 

modes  for  generating  intense  microwave  radiation.  We  further  note 

from  Fig.  8 that  the  growth  rate  of  the  TM  mode  is  comparable  with 

that  of  the  TE  mode,  when  Bg  approaches  unity. 
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Finally,  we  conclude  this  section  by  pointing  out  two  important 
areas  where  the  present  analysis  can  be  extended.  First, 
as  previously  shown  for  s=l,  it  may  also  be  possible  to  demonstrate 
analytically  for  higher  magnetic  harmonic  numbers  (s>2)  that  the  region 
of  k-space  corresponding  to  instability  is  very  narrow-band  with  kc-wB^- 
Second,  the  present  stability  analysis  can  be  extended  in  a relatively 
straightforward  manner  to  magnetic  harmonic  numbers  s>4,  and  also  to 
higher  values  of  Bq  (8q-v1). 


26 


1 

VI.  CONCLUSIONS 

In  this  paper,  we  have  examined  the  excitation  of  electromagnetic 
waveguide  modes  by  the  cyclotron  maser  instability  for  magnetic  harmonic 

number  s>l.  The  analysis  was  carried  out  for  a hollow  electron 

<"■  » 

beam  propagating  parallel  to  a uniform  magnetic  field  Bge^,  assuming 
that  the  beam  is  very  tenuous  [Eq.  (1)].  In  Sec.  II,  equilibrium 
properties  were  calculated  for  the  choice  of  electron  distribution 

function  in  which  all  electrons  have  the  same  value  of  canonical  angular 

2 

momentum  (Pg)  and  same  value  of  energy  (yqHic  ) in  a frame  of  reference 

moving  with  axial  velocity  8,c  [Eq.  (5)].  A formal  stability  analysis 

for  azimuthally  symmetric  electromagnetic  perturbation  was  carried  out 

in  Sec.  III.  Equations  (38)  and  (39),  when  combined  with  Eqs.  (16), 

(22)  , (23)  , and  (24) , constitute  one  of  the  main  results  of  this 

paper  and  can  be  used  to  investigate  stability  properties  for  a broad 

range  of  system  parameters.  In  Sec.  IV,  a detailed  analytic  investigation 

of  the  cyclotron  maser  instability  was  carried  out  for  TE  and  TM 

2 _ 2 

waveguide  modes,  assuming  3q=(1-1/Yq)<<1.  The  value  of  the  geometric 
parameter  Rg/R^  corresponding  to  maximum  growth  rate  was  determined, 
and  a detailed  numerical  analysis  of  the  dispersion  relation  was 
presented  in  Sec.  V.  One  of  the  principal  conclusions  of  this  study 
is  that  for  moderate  or  high  value  of  Eg  (Eg>0.3,  say)  magnetic 
harmonic  perturbations  with  s>2  have  growth  rates  comparable  with  t:,e 
fundamental  (s=l)  mode.  Moreover,  the  growth  rate  of  perturbations 
with  s>2  increases  rapidly  when  the  value  of  Eg  is  increased. 

Finally,  we  emphasize  that  Eq.  (53)  together  with  Eq.  (30), 


can  be  used  to  calculate  conditions  for  maximum  microwave  generation 
by  the  cyclotron  maser  instability.  By  selecting  the  value  of 
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the  applied  magnetic  field  according  to  F.q.  (30),  and  choosing  Rq/Rc= 

a , /«_  , the  TE  mode  growth  rati-  can  he  maximized  for  magnetic  harmonic 
si  On 

number  s,  thereby  optimizing  the  microwave  power  output  for  radiation 
with  frequency  ai^su^y^. 
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FIGURE  CAPTIONS 


Equilibrium  configuration  and  coordinate  system. 

Electron  density  profile  [Eq.  (13)]. 

The  straight  lines  u=kVj+sa  /y^  and  u)=kc/B^  intersect  at 

(w_  ,kri)  = (sw  y,  , sw  y,/c)  . For  the  TE  mode,  the  curve 
0 0 c b c b 

2 2 2 2 2 1/2 

u)=(k  c +a_  c /R  ) passes  through  (w„,k„)  provided 
On  c u u 

2 2 2 2 2 1/2 

an  c/R  =sw  . For  the  TM  mode,  the  curve  w=(k  c +8  c /R  ) 

On  c c Un  c 

passes  through  (c^.k^)  provided  ^qi1c'/Fc=swc* 

Electron  orbit  in  a plane  perpendicular  to  the  z-axis. 

Plots  of  (a)  normalized  TE  mode  growth  rate  y.w./u)  , 

b l c 

and  (b)  normalized  Doppler  shifted  real  frequency  Re(y^fi^/wc) 

[y,  (w  -kV,)-w  ]/u  versus  Rn/R  [Eq.  (48)]  for  s=l , v=0.001, 
b r b c c U c 

8^=0. 4,  6,  =0.286  and  several  values  of  n. 
u b 

Plots  of  normalized  TM  mode  growth  rate  v-to./w  versus  R„/R 

b i c 0 c 

[Eq.  (51)]  for  parameters  identical  to  Fig.  5. 

Plots  of  normalized  TE  mode  growth  rate  Y.w./to  versus  R„/R 

b i c 0 c 

[Eq.  (48)]  for  (a)  s=2,  and  (b)  s=3,  and  parameters 
otherwise  identical  to  Fig.  5. 

Plots  of  (a)  normalized  maximum  growth  rate  (y,w./w  ) , 

b i cm 

and  (b)  normalized  Doppler  shifted  real  frequency  Re(y^n/wc) 

versus  8_,  for  n=2,  v=0.001,  B,=8r>/1.4  and  several  values  of 
u b u 

s.  The  solid  and  broken  curves  correspond  to  the  TE  and  TM 
modes,  respectively. 





(a)  TE  mode,  i/  = O.OOI , s = I , /3=0.4  , /3h=0.286 


